Isotopic effects of hydrogen adsorption in carbon nanotubes 



M.C. Gordillo, J. Boronat, and J. Casulleras 

Departament de Fisica i Enginyeria Nuclear, Campus Nord B4-B5, 
Universitat Politecnica de Catalunya. E-08034 Barcelona, Spain 
(February 1, 2008) 



Abstract 

We present diffusion Monte Carlo calculations of D2 adsorbed inside a nar- 
row carbon nanotube. The ID D2 equation of state is reported, and the 
one-dimensional character of the adsorbed D2 is analyzed. The isotopic de- 
pendence of the constitutive properties of the quantum fluid are studied by 
comparing D2 and H2. Quantum effects due to their different masses are ob- 
served both in the energetic and the structural properties. The influence of 
the interatomic potential in one-dimensional systems is also studied by com- 
paring the properties of D2 and 4 He which have nearly the same mass but a 
sizeably different potential. The physics of molecular hydrogen adsorbed in 
the interstitial channels of a bundle of nanotubes is analyzed by means of both 
a diffusion Monte Carlo calculation and an approximate mean field method. 
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I. INTRODUCTION 



Since their discovery by Iijima in 1991,0 carbon nanotubes and their fascinating properties 
have attracted the attention of both theoretical and experimental physicistsBa Nanotubes 
are basically long cylinders, their walls being formed from one (single-walled carbon nan- 
otubes, SWCN) or several (multiple-walled ones, MWCN) graphite sheets. One of its more 
important features is its extremely narrow width, with diameters in the nanometer scale, 
ranging from ~ 7 to 40 A, compared to its length, thousands of times larger. The enormous 
length to diameter ratio and the narrowness of the tube make nanotubes excellent systems 
to hold inside nearly perfect one- dimensional fluids. Moreover, as a consequence of its nar- 
rowness, the adsorption energy of an atom or molecule inside a carbon nanotube is several 
times greater than on planar graphite. This feature ensures that any particle of adequate 
size will be swallowed inside carbon nanotubes, provided that the half-fullerene caps at their 
ends are removed. Its use as a possible solution to the packing of H 2 in fuel cells has been 
discussed in the literature.!! - ! 

Apart from its technological interest, the quasi-one dimensionality of these carbon tubes 
is itself an interesting property from a purely theoretical point of view. In the same way 
that a sheet of graphite provides an almost perfectly 2D experimental environment, the 
small transversal room available to the particles adsorbed in a narrow carbon nanotube 
creates a proper setup to extract experimental information of quasi-one-dimensional systems. 
Additionally, if the temperature is low enough, quantum effects become important, and an 
experimental realization of a quasi-one-dimensional quantum fluid is realized. A first step 
in that direction has been recently given by Teizel et a/,0 who unambiguously observed 
the quasi-one dimensional behaviour of 4 He adsorbed inside SWCN bundles. The helium 
intake could be directed either to the inside of the nanotubes, to the interstitial channels 
between the several tubes constituting a bundle, or to the external surface of the nanotube 
bundle. The occupation probability in the three regions depends on the geometry of the 
bundle, on the size of the atom/molecule to be adsorbed and, finally, on the interatomic 
potential between the carbon atom and each particular species.0 In all cases, an almost 
perfect alignment of the atoms is granted. Encouraged by this experimental accessibility, 
some theoretical work about the adsorption of quantum liquids inside nanotubes has been 
carried out. The most studied case is 4 He, which has been considered both in ID and inside 
a (5,5) nanotube. In the limit of zero temperature, both ID and quasi-one-dimensional 4 He 
are self-bound systems, with a binding energy ranging the milikelvin scale.&"0 

Molecular hydrogen inside nanotubes is also a very appealing system. In fact, the con- 
finement of H 2 in SWCN, in the limit of zero temperature, has recently deserved theoretical 
work. Its mass is approximately half the 4 He one, but the hydrogen-hydrogen interatomic 
potential has a potential well three times deeper than the helium-helium one. This prevents 
the existence of a liquid phase at zero temperature in both 2D and 3D, but diffusion Monte 
Carlo (DMC) calculations^ have shown that it is not enough to preclude it in ID. This 
result is due to the reduced number of neighbours of each atom. In fact, this liquid state 
could be produced also in higher dimensions when the number of molecules surrounding a 
particular one is artificially reduced, as in small clusters^ (3D) and in surfaces doped with 
the right kind of impurities^ (2D). 

In the present work, the DMC method is used to study the influence of both the in- 
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terparticle potential and mass on the thermodynamic behavior of the isotopes of molecular 
hydrogen, especially H 2 and D 2 , adsorbed inside a carbon nanotube. Since the electronic 
structure of molecular deuterium and hydrogen is the same, the D 2 -D 2 interparticle potential 
is identical to the H 2 -H 2 one. This is equally true for the particle-tube interactions. More- 
over, the mass of the D 2 molecule is very similar to that of a 4 He atom. Thus, by comparing 
the deuterium results with those for helium, the effect of the respective potential wells can 
be inferred. On the other hand, the influence of the zero-point energy in the thermodynamic 
behavior of a quasi-one dimensional array can be drawn from the comparison between the 
deuterium and hydrogen results. In the last part of the paper, we present results for the 
equation of state of H 2 , D 2 , and 4 He when they are adsorbed in the interstitial channels of 
a nanotube bundle. The influence of neighbouring arrays in a particular one, as a function 
of the interchannel separation, is studied using the DMC equation of state of purely ID 
systems. The attractive Van der Waals tails between the arrays increase the binding energy 
of the adsorbed liquid in a quantity which decreases with the interchannel separation. Using 
the DMC method, it is verified that a mean field approximation, as suggested in Ref. 0, is 
very accurate at realistic separations. 

The outline of the paper is as follows. In the next section, we briefly introduce the 
DMC algorithm used to study both the fluid adsorbed in a nanotube and the purely one- 
dimensional model. In Sec. Ill, we compare results of the equation of state and spatial 
structure for the two isotopes, H 2 and D 2 . Predictions for the equations of state of both 
isotopes, when the adsorbtion is in the interstitial channels, are also reported. Finally, Sec. 
IV comprises the summary and the main conclusions of our work. 



II. METHOD 

The diffusion Monte Carlo method^il is a theoretical tool that has proved its accuracy 
in a large variety of systems and physical scales. It solves in a stochastic way the Schrodinger 
equation, providing results for the ground state which are exact for boson systems like 
the present one. The iV-body Schrodinger equation is solved in imaginary time after the 
introduction of importance sampling, a standard technique to reduce the variance and make 
the algorithm operative. Explicitely, 

- = -D V 2 /(R, t) + DV(F /(R, t)) (1) 

+ (E L (R)-E)f(R,t) , 

with /(R, t) = V'(R) ^(R ; t), where ip(R) denotes a trial wave function used for importance 
sampling and \I/(R, i) the imaginary-time-dependent wave function. In Eq. (Q), D = % /2m, 
El(R) = ^(R) _1 ^(R) is the local energy, and F(R) = 2^(R)" 1 V^(R) is the so-called 
drift force which guides the diffusion movement to regions where ip(R) is large. 

DMC calculations have been carried out for two different situations: i) a pure ID array 
of D 2 molecules, and ii) D 2 filling a tube of radius R = 3.42 A, corresponding to the (5,5) 
armchair nanotube in the standard nomenclature .ill This structure has been chosen because 
it is one of the narrowest nanotubes experimentally obtained, and therefore it is expected 
to be close to an ideal ID arrangement. 
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In the present calculation, we have considered the D2 molecules interacting via the Sil- 
vera and Goldman (SG) pair potential.^ This semiempirical interaction is the same that has 
been used in Ref. [15] in the study of H 2 in tubes, and has become the standard approach in 
molecular hydrogen and molecular deuterium calculations. The SG interaction is a spher- 
ically averaged model, even though the pair of molecules involved are ellipsoids. However, 
the eccentricity of the ellipsoids is so small (r = 0.94, in the H 2 case) that it produces results 
of remarkably accuracy, even at low temperatures. For the D2-tube interaction we have used 
the Stan and Cole (SC)BJH H 2 -tube interaction. This is a cylindrically averaged potential, 
i.e., the potential felt by a molecule inside the tube is a function only of its distance to 
the center of the tube. We have verified that the differences between the results obtained 
with this averaged interaction and the ones derived by explicitely summing up all the C-D 2 
interactions inside the tube are negligible. The SC potential depends only on three param- 
eters: a, e, and the radius of the tube. For the (5,5) tube considered, R = 3.42 A, and 
the parameters a and e correspond to the Lennard- Jones C-H 2 interaction (a = 2.97 A and 
e = 42.8 K).H In Fig. 1, the H 2 -tube SC potential is displayed and compared to the 4 He 
potential inside the same cylinder. The two curves reflect the differences between the pair 
(a, e) for 4 He and H 2 ; the larger values of e and a in the case of H 2 originates a deeper well 
and a shorter range, respectively. 

The trial wave function ip used in the ID system is 

tf 1D (R) = tfj(R), (2) 

with \l/j(R) = ]lj<j ex P (—0.5 (bj r^) 5 ) a Jastrow wave function with a McMillan two-body 
correlation factor. 

For the simulation of D 2 inside the tube, an additional one-body factor is introduced in 
order to avoid the hard core of the tube-molecule interaction, 

v|> T (R) = ttj(R)* c (R) , (3) 

with \& C (R) = Y\f exp(— erf), r { being the radial distance of the particle to the center of the 
cylinder. 

The parameters b and c entering in and \l/ c , respectively, have been determined by 
means of a variational Monte Carlo (VMC) optimization. Near the equilibrium density 
we find the optimal value b = 3.996 A, gradually increasing with the linear density A (at 
A = 0.358 A -1 , b = 4.026 A). The parameter c, inversely proportional to the width of the 
gaussian defined in \l/ c , has been fixed to c = 6.392 A -2 due to its negligible dependence on 
the density. It is worth mentioning that these parameters are different from the ones used 
for H 2 in the same setupEl (b = 3.759 A, c = 4.908 A" 2 ). 

III. RESULTS 

A. Low-density regime 

In order to study the behavior of the several isotopes of molecular hydrogen inside nan- 
otubes, the dependence of the adsorption energy on the mass species constitutes a first 
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relevant point. Since the interparticle and tube-molecule potentials do not distinguish be- 
tween isotopes, their different masses play the central role. At low densities, the dominant 
effect is due to the binding energy of single independent molecules. From our DMC calcu- 
lations the binding energy of a single molecule inside a carbon tube turns out to be, for the 
(5,5) tube, -1539.87 ± 0.11 K for H 2 , and -1605.23 ± 0.01 K and -1624.37 ± 0.01 K for 
D2 and T2, respectively. The increase in the binding energy with the mass comes from the 
combination of two features: a decrease in the kinetic energy, mainly due to a direct effect of 
the mass (m^/m^ ~ 1/2, m-o 2 /m^ 2 ~ 2/3), and a simultaneous increase of the potential 
energy. Although those energies correspond to the ground state of single molecules at K, 
they constitute a very good estimation in the limit of infinite dilution at nonzero tempera- 
tures. From the above binding energies one can extract information on the selectivity in the 
adsorption inside the nanotube. Following Ref. the selectivity of isotope 2 with respect 
to isotope 1 can be defined by the quotient S = {xx/x^)/{yx/y^) with X, (y^) the nanotube 
(bulk) mole fractions. It has been proved that in the limit of zero pressure the selectivity 
So is very well approximated by 

m 2 / E x -E 2 \ 
S = — exp , (4) 



mi V T 

where E^ is the binding energy of isotope i. Considering T = 20 K, as in Ref. |24|, we obtain 
S , o(T 2 /H 2 ) = 22.8 and S , (T 2 /D 2 ) = 1.7 for the (5,5) tube. The selectivity is especially high 
in the case T 2 /H 2 due to the sizeable difference in binding energies between the two isotopes, 
Et 2 — E\\ 2 = —84.5 K. That large selectivity, which is a purely quantum effect, has been 



proposed in Ref. |23] to achieve an efficient isotopic sieving. 



B. One- and quasi-one-dimensional systems 

The ground-state properties of an array of hydrogen molecules, in a one- and in a quasi- 
one-dimensional environments, have also been studied using the DMC method. In a previous 
work,E3 we had studied the case of H 2 : the equilibrium state in ID corresponds to a liquid 
phase with a density A = 0.2191 ± 0.0004 A^ 1 and an energy per particle (E/N)q = 
—4.834 ± 0.007 K. A first approach to determine the equation of state of ID D 2 is based 
on the use of the average correlation approximation (ACA). This approximation, which has 
been widely used in the 3 He- 4 He isotopic mixture in the limit of zero 3 He concentration,!!!! 
consists in the present case in evaluating the expectation value of the Hamiltonian for D 2 
using the exact wave function of the other isotope (H 2 ). Thus, the energy per particle of D 2 
in ACA is given by 



D 2 




(5) 



where (T/N)n 2 is the kinetic energy of the H 2 molecule. 

ACA corresponds, obviously, to a variational estimate and thus the D 2 energies obtained 
from Eq. (M) are upper bounds to the exact values. Using the DMC kinetic energies of ID 



H 2 reported in Table I, and the total energies published in Ref. 15, we have calculated the 



ACA equation of state for D 2 . The results obtained are shown in Fig. 2 and compared with 
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the exact solution derived from a DMC calculation of ID D2. Both equations of state (lines) 
correspond to a least squares fit to the data in the usual form 



The parameters (E/N) , X , A, and B are given in Table II. The upper bounds provided 
by ACA appear clearly reflected in Fig. 2, its quality being quite reasonable taking into 
account its straightforward estimation (|5]). In fact, the discrepancies between ACA and the 
exact results are less than 2 % for the equilibrium density value, and around 3 % for the 
binding energies. ACA can also be used to estimate the equation of state of the T2 fluid, a 
system that we have not studied here using the explicit DMC calculation. For ID T2, ACA 
predicts A = 0.2501 A" 1 and (E/N) = -12.677 K. 

On the other hand, the equation of state of D 2 inside a tube can also be estimated by 
ACA, since the tube-molecule SC potential does not depend on the particular isotope. Using 
the H2 equation of state (Ref. [15]) and the H2 kinetic energies reported in Table I, the ACA 
results compare with the corresponding DMC ones with an accuracy similar to the ID case. 
DMC energy results for both ID D 2 and D 2 adsorbed in the (5,5) nanotube are displayed 
in Fig. 3. In order to show the two equations of state with the same energy scale, we have 
subtracted to the tube results the adsorption energy of a single molecule. In Fig. 3, the 
curves are polynomial fits (^) to the DMC data, the optimal parameters for the tube being 
reported in Table II. It can be seen that the equilibrium density for ID D2 and D2 adsorbed 
in the (5,5) nanotube are almost identical. This is also true for the location of the spinodal 
points of D 2 , A S 1D = 0.230 ± 0.001 A" 1 and Xj = 0.232 ± 0.001 A -1 , which can be derived 
from the data contained in Table II. 

The binding energies, in the respective equilibrium points ((E/N)o), are slightly different: 
the additional transverse degree of freedom only amounts to an increase of 0.091 K. This 
increase in the binding energy is nearly a factor two smaller than the one drawn from the 
DMC calculations for H 2 (0.172 K). In relative terms, the increase of the binding energy is 
only a 0.85 % for D 2 versus a 3.5 % for H 2 . Therefore, one can conclude that the effects 
of the additional degree of freedom of the D2 molecules in the radial direction inside the 
nanotube, which account for the enhancement of the binding energy, are reduced by the 
greater mass of the D2 molecule with respect to the H2 one. As a matter of comparison, 
it is illustrative to compare the effects observed in D 2 with the ones previously studied in 
4 He using the same methodology and geometry. It is worth noticing that the masses of 
D 2 and 4 He are nearly the same whereas the interatomic potentials are sizeably different. 
The DMC results show that the latter effect is completely determinant: in 4 He the relative 
difference mentioned above is 90 %, two orders of magnitude larger than in D2. Another 
minor effect that contributes to the one- dimensionality of molecular deuterium adsorbed 
inside the tube, is the larger hard-core size of the C-D 2 interaction (o"c-d 2 — 2.97 A), versus 
the C-He one (<Tc-He — 2.74 A). The mass versus potential effects can also be seen in the 
value of the equilibrium density. Inside the tube, Ao goes from 0.079 ± 0.003 A" 1 in 4 He, 
to 0.2200 ± 0.0006 A" 1 in H 2 , to reach 0.2473 ± 0.0002 A" 1 in D 2 . That sequence clearly 
indicates that the main influence in Ao comes from the potential, since the isotopic change 
varies the location of the energy minimum less than 15 %. Most probably, the features 
observed in these systems have a more general character and one can guess that if the well 
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of the interatomic potential is increased and/or the mass of the particle adsorbed inside a 
tube is enlarged, the effect would be an increase in the Ao value. 

In Fig. 4, the density dependence of the pressure for both ID D 2 and D 2 inside the tube is 
reported from the equilibrium density up to 0.30 A -1 . For the sake of comparison, the same 
results for H 2 are also plotted. In both H2 and D 2 , the pressures for the ID systems (solid 
lines) are greater than the ones for the tubes (dashed lines) at the same A, the difference 
being quite similar for the two isotopes. For the same geometry and at equal density, the 
pressure of H 2 is always larger, a fact that is basically due to its slightly smaller equilibrium 
density. 

Complementary and useful information on the system can be obtained from the micro- 
scopic study of the spatial structure of the molecules in the array. In Fig. 5 the radial 
distribution functions, g z {f), along the z axis, are shown. They correspond to the quantum 
fluids adsorbed in the tube at their respective equilibrium densities A . The ID counterparts 
are nearly identical and are not displayed. Being the denser of the three systems, D 2 exhibits 
accordingly the most pronounced oscillations in the g z (r) function. The shift in the positions 
of the maxima for the two molecular isotopes arises basically from the difference in their 
respective Ao's. A comparison of the radial distributions of H 2 and D 2 at the same linear 
density, not shown for simplicity, indicates so. In the 4 He case, the much lower equilibrium 
density, direct consequence of the different potential, explains the smoothness of the g z {r) 
obtained. 

The radial densities inside the (5,5) tube have been also studied. In Fig. 6, the radial 
densities for 4 He, H 2 , and D 2 for the same linear density A = 0.245 A -1 are shown. The 
trends shown in the figure are common to all densities studied: the particle with the largest 
mass (D 2 ) is the one which spends more time in regions closer to the center of the tube, i.e., 
D 2 in the tube is the closest to a one-dimensional system. The change in the mass and in 
the interatomic potential work in the same direction: the radial densities of H 2 and 4 He are 
quite similar. Both curves show a decrease in the radial localization and larger fluctuations 
in the transverse direction. 



C. Adsorbtion in a bundle of nanotubes 

The absorption of different species inside the interstitial channels of bundles, formed when 
several nanotubes are staked together, has been theoretically predicted^ and experimentally 
observed.0 The staking of the nanotubes allows for the formation of an hexagonal network 
of quasi-one dimensional channels separated by distances which depend on the diameter 
of the nanotubes themselves. For instance, when the nanotube ropes are made out of 
(5,5) constituents, the distance between two adjacent tubes is 10.26 A, whereas the closest 
interstitial channels are located ~ 5.9 A apart. The interstitial channels are very narrow, 
with effective diameters oscillating between 3 and 3.5 A.H This size is considerably smaller 
than the diameter (6.8 A) of the (5,5) tube considered in the previous subsection. Therefore, 
it is plausible to assume that if molecular hydrogen or deuterium is introduced inside one of 
this channels, its behavior would have even a more one dimensional character than in the 
inner part of a single tube. 

Considering the previous arguments, a quantum fluid adsorbed in the intersites can be 
well modelled by hexagonal arrangements of purely ID systems, mutually separated by the 
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real interchannel distance in the bundle. To a large extent, the additional consideration of 
the interactions between the adsorbed atoms/molecules and the carbon atoms of the nan- 
otube walls would simply introduce a constant shift in the energy scale. The influence of 
neighbouring channels in the binding energy of 4 He adsorbed in one of them has been esti- 



mated in Ref. [18| considering a mean field approximation. In this approach, the neighbouring 
channels are seen as uniform arrays that only contribute to increase the potential energy in 
the form 

\ roc , . 

A m{ = -j_^dxV(V^Td?) , (7) 

where d is the distance between the reference channel and one of its neighbours. The total 
energy correction is obtained by summing up the contribution of the nearest neighbours, 
the next-nearest neighbours, and then on up to a desired accuracy. From its definition (0) 
it is clear that the energy per particle in this approximation, 

E/N= (E/N) 1D + A mi , (8) 

is an upper bound to the exact energy. 

In order to test the mean field approximation (^), we have performed full DMC calcu- 
lations of a system formed by an hexagonal array of ID H 2 channels. Since the distances 
between channels depend on the diameter of the tubes, the calculations have been made in 
the more challenging case for Eq. @, a bundle of (5,5) nanotubes. That tube is one of 
the narrowest experimentally obtained, and therefore the interchannel distance is one of the 
smallest in nature. Actually, that is an unrealistic case since the corresponding interstitial 
channels are too narrow to adsorb H 2 . The DMC energies are displayed in Fig. 7, and com- 
pared with the approximated estimations @. As it is evident in the figure, both results are 
nearly identical, at least in the density range considered. The conclusion is that correlations 
between hydrogen molecules located in different channels suppose a negligible contribution 
to the energy, even in this case where the channels are very close. In Fig. 7, there is also 
reported the equation of state of ID H 2 . As one can see, the increase of the binding energy 
of the system, when the array is immersed in a bundle, is as large as a factor two. 

Having checked the validity of the mean field approach (§) in an exigent case, one can 
use that approximation to study the equation of state of molecular hydrogen and helium 
adsorbed in other more appropriate bundles. It is worth mentioning that the number of 
neighbours to be included in A m f has to be larger than ten, otherwise the binding energy is 
underestimated in a 5 %. Besides the (5,5) bundle, with an interchannel distance a = 5.92 A, 
we have analyzed the following cases: (6,6), a = 6.71 A; (8,8), a = 8.27 A; (10,10), a = 9.83 
A. The fluids that have been analyzed are 4 He, D 2 , and H 2 . In Fig. 8, the differences 
between the equilibrium density in the bundles and in the purely ID system are shown for 
the three cases. There is a clear difference (an order of magnitude) between the increase of 
Ao, when 4 He is adsorbed in the intersites, and the density shift experimented by the two 
hydrogen isotopes. In H 2 and D 2 the change in density is small but one can distinguish 
a slightly larger effect in the lighter isotope. In all cases, the density difference decreases 
when the interchannel separation increases. However, H 2 and D 2 have already reached the 
ID equilibrium density in the (10,10) geometry whereas 4 He still shows a sizeable effect. 

The increase of the binding energies with respect to the one-dimensional systems are 
reported in Fig. 9. The different points correspond to the same bundles reported in Fig. 8. 
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In absolute terms, that increase is much larger in molecular hydrogen than in 4 He, the effect 
being slightly larger in D 2 than in H 2 for all the bundles. The energy difference between H 2 
and D 2 is very small and decreases when the interchannel distance increases. It is clear, from 
the mean field correction (0), that the difference in energy between the two isotopes arises 
from their different densities since the interatomic potential is the same. In 4 He the shift 
in energy is comparatively small but, taking into account that the ID 4 He binding energy^! 
is only (E/N)q = —0.0036 K, one concludes that the bundle effects are much larger in 4 He 
than in molecular hydrogen. 



IV. SUMMARY AND CONCLUSIONS 

In this work, we have studied the zero-temperature equation of state of molecular deu- 
terium, and compared the results with those previously obtainedEl for 4 He and molecular 
hydrogenS The results show that the different masses of the two isotopes do influence both 
the energetic and structural properties. As a general trend, D 2 appears as a slightly denser 
liquid, a fact that is observed in its equation of state and also in the radial distribution 
functions. The comparison between D 2 and 4 He, two particles with the same mass but with 
very different interatomic potentials, has shed light on the influence of the potential in one- 
dimensional systems. Concerning the mass and potential changes, the DMC results show 
unambiguously that the potential effects are manifestly dominant. Bundles of carbon nan- 
otubes have also been studied by modelling them as hexagonal one-dimensional arrays. The 
explicit DMC calculation of the system has allowed to check the accuracy of a previously 
suggested mean field approximation]^ Once verified the validity of that approximation, we 
have analyzed the effects of the neighbouring channels on the ID equation of state. Our 
results show that those effects are much more relevant in He than in molecular hydrogen. 

In the tube and bundle calculations, we have assumed that the interparticle potential 
is not modified by the existence of nearby carbon atoms. Recently, Cole and coworkersilil 
have estimated the effective correction in the two-body potential due to three-body inter- 
actions between two particles and the carbon atoms of the nanotube wall. This three-body 
interaction has been modelled by the well-known triple-dipole interaction derived by Ax- 
ilrod and Teller,!! a potential energy contribution that is mostly repulsive. In the case of 
molecular hydrogen in the interstitial channels, the well of the effective Silvera and Goldman 
potential is reduced by approximately a factor two. This unexpected and big effect would 
modify some of the results presented in our work. Nevertheless, it is well-known that, in 
spite of the fact that Axilrod- Teller is the dominant contribution to V3, at short interparticle 
distances an attractive force emerges. This short-ranged three-body interaction, known as 
exchange interaction]^ is due to the influence of a third particle in the charge densities of 
two interacting atoms. Due to the short interparticle distances H 2 -C-H 2 in the channel, that 
attractive term would partially cancel the repulsion introduced by Axilrod- Teller. We have 
preferred to ignore the three-body corrections!! until a better knowledge of the nonadditive 
terms of this particular system is achieved. 

We plan, in the near future, to introduce in the simulation the possibility of the swelling 
of the channels when molecular hydrogen is adsorbed. The aim would be to test, by a 
microscopic calculation, the recent theoretical result!! that H 2 forces the nanotube bundles 
to slightly spread upon adsorption. We hope our calculations, and the ones carried out by 
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other theoretical groups, can stimulate new experiments on molecular hydrogen adsorption 
in carbon nanotubes and bundles, especially at low-temperature where quantum effects 
become macroscopic. 
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TABLES 







ID 








Tube 




A (A-i) 


H 2 




D 2 


H 2 




D 2 




0.225 


9.36 ± 0.03 




5.63 ± 0.02 


125.5 ± 


0.2 


80.5 ± 


0.1 


0.230 


10.20 ± 0.03 




6.03 ± 0.02 


125.8 ± 


0.2 


80.9 ± 


0.1 


0.235 


11.32 ± 0.02 




6.57 ± 0.02 


127.1 ± 


0.2 


81.2 ± 


0.1 


0.239 


12.30 ± 0.03 




7.28 ± 0.02 


127.8 ± 


0.2 


81.8 ± 


0.1 


0.259 


18.27 ± 0.05 




10.98 ± 0.03 


133.6 ± 


0.2 


85.4 ± 


0.1 


TABLE I. Kinetic energies 


of H 2 


and D 2 in ID and inside 


: a (5,5) carb 


on nanotube. 







ACA 


ID D 2 


D2 in a tube 


Ao (A-i) 


0.2419 ± 0.0007 


0.2457 ± 0.0003 


0.2473 ± 0.0002 


(E/N) (K) 


-10.282 ± 0.008 


-10.622 ± 0.016 


-1615.94 ± 0.015 


A(K) 


154 ± 3 


2.0 10 2 ± 1.0 10 1 


2.13 10 2 ± 1.0 10 1 


B (K) 


4.5 10 2 ± 70 


9.6 10 2 ± 1.2 10 2 


1.10 10 3 ± 1.1 10 2 



TABLE II. Parameters of the equation of state of D 2 (Eq. 
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FIGURES 




FIG. 1. SC tube- molecule/atom potentials inside a (5,5) carbon nanotube. 




-10.8 b 1 1 1 1 1 1 1 

0.23 0.24 0.25 0.26 

X (A" 1 ) 

FIG. 2. Equation of state of ID D2. The symbols are the DMC energies and the solid line a 
polynomial fit to them (||). The dashed line corresponds to the energies estimated using ACA. 
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FIG. 3. Comparison between the equation of state of ID D2 and D2 adsorbed in the nanotube. 
Filled squares correspond to the tube results; open squares, to the ID ones. The lines are poly- 
nomial fits to the DMC data. To better compare both results, we have subtracted to the tube 
energies the binding energy of a single molecule. 




FIG. 4. Pressure of D2 and H2 as a function of the linear density A. Solid and dashed lines cor- 
respond to the ID and nanotube systems, respectively. The left (right) scale is the ID (nanotube) 
pressure. 
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FIG. 5. Two-body distribution function in the nanotube system, and in the z direction. Solid, 
dashed, and dotted lines correspond to D2, H2, and 4 He, respectively. All curves are calculated at 
their respective equilibrium densities Aq. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

r(A) 

FIG. 6. Radial density of D2 (solid line), H2 (dashed line), and 4 He (dotted line) inside the 
(5,5) nanotube. 
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FIG. 7. Equation of state of D2 in a bundle of (5,5) nanotubes. The symbols are the DMC 
results, and the dotted line corresponds to the mean field approximation. The dashed line is the 
ID result. 
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FIG. 8. Difference between the equilibrium densities, in a bundle and in a ID array, as a 
function of the cell parameter a. Diamonds, filled squares, and open squares correspond to D2, H2, 
and 4 He, respectively. 
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FIG. 9. Difference between the binding energies, in a bundle and in a ID array, as a function of 
the cell parameter a. Diamonds, filled squares, and open squares correspond to D2, H2, and 4 He, 
respectively. 
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